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Polygonization of carbon nanotubes 
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We use a multiscale procedure to derive a simple continuum model of multiwalled carbon nan- 
otubes that takes into account both strong covalent bonds within graphene layers and weak bonds 
between atoms in different layers. The model predicts polygonization of crossections of large mul- 
tiwalled nanotubes as a consequence of their curvature-induced turbostratic structure. 
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Polygonization of crosssections of large multiwalled 
carbon nanotubes (MWNTs), especially following heat 
treatment, has been observed in a number of studies 
0) S HI- The physics of this phenomenon is gener- 
ally well-understood — the curvature-induced mismatch 
between the lattices of the adjacent graphene shells can 
be mollified by flattening the shells at the expense of cre- 
ating the line defects needed to maintain the cylindrical 
shape of the tube 

The modeling of polygonization have proved more elu- 
sive: the atomistic modeling of large tubes is computa- 
tionally expensive, while continuum models that do not 
take into account atomic structure of graphene are inade- 
uate to reproduce lattice-dependent features of MWNTs 
To our knowledge, the only existing model Q 
of polygonization in MWNTs is rather complex and it 
assumes from the start that a nanotube is polygonal. 
The extensive first-principle calculations are then used 
to compute the energy of the defects and the planar sec- 
tions of the tube in order to select the optimal shape of 
MWNT. 

In this letter we employ a continuum theory derived 
by upscaling a simple atomistic model to predict poly- 
gonization of large MWNTs. Our main premise is that 
polygonization is the result of a competition between the 
two types of interatomic interactions — the interactions 
within individual graphene layers and the interactions 
between these layers — and it is controlled by the diame- 
ter but not the number of walls in the MWNT. We take 
advantage of the last observation by considering the sim- 
plest case of a two- walled carbon nanotube. 

To derive a continuum model of a two- walled nanotube 
we will assume that the diameter of each tube is much 
larger than the length of an interatomic bond in a hexag- 
onal carbon lattice. Then a microscale object has dimen- 
sions of order of the carbon bond length and a macroscale 
object has dimensions of order of the diameter of the nan- 
otube. 

Macroscopic geometry: Suppose that two smooth, 
closed concentric curves C\ and C 2 represent a crossection 
of a two-walled carbon nanotube by a plane perpendicu- 
lar to the axis of the tube (Fig. [1]). To simplify mathe- 
matics we will assume that these curves are parallel, that 



is the distance between C\ and C 2 when measured along 
any normal to C\ is equal to the same constant value of 
d. The R 2 -valued function fi(s) parametrizes C\ with 
respect to its arclength s 6 [0,L] measured counterclock- 
wise from some fixed initial point r(0). Here L is the 
length of Ci . 




FIG. 1: Crossection of a two- walled carbon nanotube. 

An orthogonal frame at f i (s) is given by a pair 
(T(s),N(s)) with T(s) = -r[{§) = (cos 0(s), sin 0(s)) 
and N(s) = (— sin0(s), cos 0(5)), where 0(5) is the angle 
between T(s) and the a;— axis at the point f i (s). The cur- 
vature of C\ can be computed via the relation k(5) = 0'(s) 
then, according to the Frenet formulas, T'(S) = k(s)N(s) 
and N'(s) = -k(s)T(s). 

In what follows, we will require that dn(s) < 1 for 
all s 6 [0, L\. This assumption is physically reasonable 
as it automatically holds for any nanotube with a con- 
vex crossection and it only rules out relatively large in- 
ward folds of the curve C\. Then we can parametrize 
C 2 with respect to s by setting f 2 (s) = ?i(s) + dN(s). 
Note that s is not an arclength parameter for f 2 . Indeed, 
fi(^) — — (1 + dre(s)) T(s) and the distance traveled 
along C 2 from ?2(0) to r 2 (s) is l(s) = f*(l + dn(a))da = 
s + (0(3) — 0(O))d while the distance traveled along C\ 
from ?i(0) to ?i(s) is equal to s. Then the running dif- 
ference between the distances traveled along C 2 and C\ is 
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C(s) — l(s) —s = d(8(s) — 0(0)) and the overall difference 
between the lengths of the two curves is C(L) = 2ird. 
Further, by an appropriate rotation of coordinates, we 
can set 0(0) = to obtain 



C(s) = dd(s) 



(1) 



Microscopic structure: Having established the macro- 
scopic geometrical framework for our model we need to 
endow it with the atomic structure of carbon. Since 
for mathematical simplicity we operate with crossec- 
tions of nanotubes that are intrinsically one-dimensional 
while the hexagonal atomic lattice of graphene is two- 
dimensional, we will assume that the projection of the 
lattice on each crossection is a one-dimensional chain of 
atoms that retains the important physical characteristics 
of the original lattice. Although this assumption might 
not lead to a quantitatively accurate effective continuum 
model, we expect that the effective model based on phys- 
ically reasonable assumptions should be qualitatively ac- 
curate. 

Suppose that two sets (chains) of equidistant carbon 
atoms are imbedded in the curves C\ and C2 , respectively, 
as shown in Fig. [2j The distance between the neighboring 
atoms in each chain will be fixed to the same constant h 
to reflect the fact that the atoms within a graphene layer 
are connected by essentially inextensible strong sp 2 a- 
bonds. Then the curve C2 contains n 
that the curve C\ . 
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Further, let the weak interaction between the atoms in C\ 
and C2 be described by a Lennard-Jones-type potential 
s( r ) = (r) 12 — (r) 6 - Then the total energy of the system 
is 
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where a and 7 are the dimensional scaling factors and N 
is the number of atoms imbedded in C\. 

From now on we will assume that d,h L and n — 
^j^- > 2 is a fixed positive integer. Nondimesionalizing 
the variables r = r/L, s — s/L and rescaling E by 7, the 
nondimensional energy is 
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+ a 



+ 



N+n-l 

E 

JV-l N+n-l 

E E 

i=i j=i 



/(tt 



r i, 



(2) 



r 2, 



where 8 = d/L « 1, 6 = h/L ^ 1, and a = a/7 are the 
nondimensional parameters of the system. Further, we 
will set Ci and C2 to represent C\ and C2 in nondimen- 
sional variables. 

Effective model: Using the smoothness of ri we have 
that — B\,i = 6'(si)e + o (e) for every i = 1, . . . , N, 

where Sj denotes the position of the z-th atom on the 
curve C\. Since / is smooth, the symmetry of / and 
the fact that it has a minimum at 8 = tt imply that 
/(tt + m -0 m+1 ) = ±/»(0'( Sl ))V + O (e 2 ) for every 
i = l,...,N. Then 
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= ^£(^ + 0(6). 



(3) 



FIG. 2: Positions of atoms on the curves Ci and C2. 

Energy: We will assume that the overall energy of the 
two-curve-system consists of two parts: the energy due 
to bending of the adjacent bonds and the energy due to 
weak delocalized 7r-bonds between the atoms imbedded 
in Ci and C2, respectively. 

Suppose that the energy associated with two bonds 
joined at an angle 6 is f(0) where / : (0, 2ir) — > R is a 
smooth convex nondimensional function that has a mini- 
mum value of at 9 = tt and satisfies f(ir — 6) = f(ir + 0). 



Since exactly the same argument applies to the analogous 
sum over the curve C2, we have that the leading contri- 
bution to the total energy due to bending of interatomic 
bonds is 



E b = fi€ (O'fds, 



(4) 



where /1 = af"(ir). This expression (cf. e.g. 0] ) corre- 
sponds to the well-known Euler elastica model for long 
slender beams @. 

Next we develop the effective expression for the weak 
interaction energy. Fix an atom i on an inner curve C\ 
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and let A — O (e 1 ^ 2 ) be an arbitrary number that we can 
use as a "mesoscale" because t « i < 1 . Set to 
be the set of all atoms on the curve C2 that are closer 
than A to the atom i. Then ^=1^9 f ^ 1 '*/ 2 '' 1 ) = 

E jeA . 9 ( J —^-) +E^a a 9 (^tH ■ The second 
sum can be approximated as follows 

h^a v / \ / 

because 

6/e = d/h = n/2n. (6) 

Having estimated the contribution to the weak inter- 
action energy from the atoms on C 2 that are macroscop- 
ically distant from the atom i 6 C\ , we now estimate the 
contribution to the energy due to the atoms that are mi- 
croscopically close to the atom i. Since the curve C 2 is 
smooth and since A <C 1, the part of the curve C2 that 
is closer than A to the atom i can be approximated by 
the tangent line to C2 at the point Sj. That is, locally, we 
have a situation depicted in Fig. [3] — in the macroscop- 
ically small neighborhood of the atom i the curves ap- 
pear as two infinite straight lines of atoms offset by some 
distance k when viewed from the microscopic perspec- 
tive. Therefore the energy of weak interaction between 
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FIG. 3: Two offset chains of atoms, 
the atom i and the atoms in is 

jeA A v ' j=-oo \ J 

to the leading order in e. 

Note that, up to a term of order e, the offset k between 
the two lattices at the point Si is equal to the running 
difference £(sj) between the two arclengths as measured 
from the points i"i(0) and r2(0), respectively. Using the 
equations ([TJ) and ^ we obtain that 



where 




Finally, 

Ei=i Ej=i g ( s - )= G(«(ai),n) 
+o(l) = \ Jo G(0(a), n)ds + o(l) = ^ + o(l), 

and the leading contribution to the total effective energy 
is given by a Ginzburg-Landau-type expression 

E eff = E b + E w = J (efi{0'f + ^G(6,n)^ds, (8) 

where e can be viewed as an analog of the Ginzburg- 
Landau parameter. 

Next we discuss the properties of the function G. Al- 
though G is defined for all of 9 and n, only integer values 
of n are physically relevant as n is equal to the difference 
between the number of atoms on the outer curve C2 and 
the inner curve C\. Then, by construction, the function 
G is 27r/n-periodic in 9 and it has exactly n minima on 
the interval [0, 2ir\. The surface of G and its crossections 
for various values of n are shown in Figs. [HO 
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FIG. 4: Surface of the effective potential G(6,n). 
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FIG. 5: G(9, n) for the different values of n. 

To summarize, in our effective model a crossection of a 
two- walled nanotube is represented by a single curve C = 
C\ with the effective energy given by 1(5)). The equilibrium 
shape of the nanotube then minimizes the energy E e ff 
subject to the constraints 

[ sin9ds= f cos6>ds = 0, (9) 
Jo Jo 

on the function 9 that enforce the closedness of C. 
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It is easy to see that the values of the energy minimiz- 
ing function 8 will reside mostly at the minima of G(9, n) 
as dictated by the "penalty" term -G(6, n). The func- 
tion 9 transitions between its two constant values over 
the narrow " interfacial" region. The width of this region 
decreases with the Ginzburg-Landau parameter e due to 
the factor of e in front of the gradient term in the energy 
E e f f . Because 9 increments by 2tt as the curve C is tra- 
versed in the counterclockwise direction and because the 
potential function has n minima on the interval [0,27r], 
the minimizing configuration of 9 should have exactly n 
interfaces separating n regions where 9 is almost con- 
stant. Note that n is equal to the number of "extra" 
atoms on the circumference of the outer tube and it is 
independent of the diameter of the inner tube (as long 
as the difference between the diameters of two tubes re- 
mains constant). 

Since e — s/L we have that e — > when L — > oo then 
the transition between the regions of constant 6 must be 
sharper for the tubes with the larger radii. On the other 
hand, the interfacial regions smear out when L is small 
and the gradient term in the energy dominates. This ef- 
fect has a clear physical explanation — the energy due to 
the curvature-induced mismatch between the lattices of 
the inner and the outer tubes increases linearly with the 
diameter while the energy of a dislocation associated with 
a corner is independent of the diameter of the tube. The 
number of dislocations needed to remove the curvature 
and to introduce the appropriate lattice stacking is equal 
to the difference between the number of atoms along cir- 
cumferences of the outer and the inner tubes. (Each dis- 
location " resets" the lattice to incorporate an extra atom 
into the outer tube.) We conclude that there must be a 
critical diameter above which polygonal crossections of 
the tubes should be energetically preferred. 

If we assume that the minimizer has the symmetry 
9(tt + s) = 2tt — 9(s), then the constraints (|9|) are auto- 
matically satisfied and the minimizer solves the boundary 
value problem 

f 2^9" - G e {6, n) = if < s < tt, , - 
\ 9(0) = and 9(tt) = tt. [ > 

This problem was solved numerically using the MATLAB 
BVP solver [3] and assuming that C2 contains seven more 
atoms than C\. The shapes of nanotubes of different 
circumferences are shown in Fig. [6] and correspond to 
the trends observed experimentally [lj. 

Observe that in our model there is no sharp transition 
between the polygonal and the circular nanotube shapes 
and the polygonal crossections continuously morph into 
the circular ones as the radii of the tubes are decreased. 
This is the consequence of the fact that the interfaces in 
the "diffuse model" are not sharp — rather, the width of 
the interfacial regions increases continuously with e. 

A simple geometrical calculation shows that, if two 
concentric circles are replaced by two concentric polygons 
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FIG. 6: Shape of a nanotube as a function of its diameter. 

of the same circumference, then the distance between the 
circles is larger than the distance between the polygons. 
However, the experimental data [H shows that the inter- 
shell spacing c/002 decreases from 0.39 to 0.34 as the diam- 
eter of the tube increases. This is the clear consequence 
of turbostraticity: for large tube diameters, polygoniza- 
tion reduces turbostraticity for the most of the nanotube 
circumference and the intershell spacing decreases to the 
equilibrium value of the correlated stacking (the spac- 
ing can be reduced further to that of the perfect AB- 
stacking if a nanotube is subjected to a high-temperature 
heat treatment .3']). The variation of intershell spacings 
for the entire range of nanotube diameters is of order 
0.39 — 0.34 = 0.05 nm — smaller than the carbon bond 
length of 0.142 nm. As a consequence, the integer part 
of the ratio of d/h remains the same for all diameters, 
validating the assumption that n (or d) can be held fixed 
in the minimization procedure. 

The extension of our approach to MWNTs in three di- 
mensions is relatively straightforward but more technical; 
it will be addressed a forthcoming paper. 
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* Electronic address: dmitry@math.uakron.edu 
[1] C.-H. Kiang, M. Endo, P. M. Ajayan, G. Dresselhaus, and 

M. S. Dresselhaus, Phys. Rev. Lett. 81, 1869 (1998). 
[2] F. Y. Wu and H. M. Cheng, J. Phys. D 38, 4302 (2005). 
[3] M. Yoon, J. Howe, G. Tibbetts, G. Eres, and Z. Zhang, 

Phys. Rev. B 75, 165402 (2007). 
[4] B. I. Yakobson, C. J. Brabec, and J. Bernholc, Phys. Rev. 

Lett. 76, 2511 (1996). 
[5] M. Arroyo and T. Belytschko, Phys. Rev. B 69, 115415 

(2004). 

[6] S. S. Antman, Nonlinear problems of elasticity, vol. 107 
of Applied Mathematical Sciences (Springer, New York, 
2005), 2nd ed. 

[7] MATLAB is a registered trademark of The MathWorks, 
Inc. http://www.mathworks.com. 



